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The band structure of cubic inverse perovskites, CasPbO and its family, are investigated with the first- 
principles method. A close observation of the band structure reveals that six equivalent Dirac electrons with 
a very small mass exist on the line connecting the P- and X-points, and at the symmetrically equivalent 
points in the Brillouin zone. The discovered Dirac electrons are three-dimensional and remarkably located 
exactly at the Fermi energy. A tight-binding model describing the low-energy band structure is also con- 
structed and used to discuss the origin of the Dirac electrons in this material. Materials related to CasPbO 
are also studied, and some design principles for the Dirac electrons in this series of materials are proposed. 
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"Emergence" is one of the most important concepts 
in condensed matter physics. Although this key word 
often appears in the context of many-body or strong 
correlation effects, emergent behaviors are also observed 
in noninteracting systems where the low-energy effec- 
tive Hamiltonian becomes quite distinct from the orig- 
inal Hamiltonian. The most well-known example is the 
relativistic Dirac Hamiltonian realized in graphene 1 ' de- 
rived from a nonrelativistic Hamiltonian. Many intrigu- 
ing properties of graphene can be ascribed to the exis- 
tence of "Dirac electrons" in its low-energy band struc- 
ture. 2 ) Actually, Dirac electrons in materials have a 
long history starting from bismuth, which has three- 
dimensional massive Dirac electrons in its band struc- 
ture. 3 ' The organic conductor a-(BETT-TTF) 2 I 3 is also 
known to be a material having Dirac electrons near the 
Fermi energy. 4 ' 5 ' The most up-to-date example is a sur- 
face state of a three-dimensional topological insulator, 6 ' 
which is extensively studied in these days. 

In connection with topological insulators, inverse- 
perovskite materials have attracted much attention re- 
cently. For example, it is claimed that CasNBi enters 
a topological phase under an appropriate strain engi- 
neering scheme. 7 ' In this paper, we show that cubic in- 
verse perovskites, CaaPbO and its family, have three- 
dimensional Dirac electrons with a very small mass at 
the Fermi energy. Although CaaPbO is on the list of po- 
tential topological insulators proposed by Klintenberg, 8 ' 
our close observation of its band structure reveals the ex- 
istence of bulk (not surface) Dirac electrons on the line 
connecting the T- and X-points, and at the symmetrically 
equivalent points in the Brillouin zone. Although some 
first-principles calculations on this material are available 
in the literature, 9 ' 10 ' it is the first time that the existence 
of Dirac electrons is pointed out. We also construct a 
tight-binding model that describes the low-energy band 
structure, and clarify the origin of the Dirac electrons in 
this material by analyzing the model. We also study the 
family of CaaPbO, and some design principles for Dirac 
electrons in this series of materials are proposed on the 




Fig. 1. (Color online) (a) Cubic inverse-perovskite structure of 
CaaPbO. (b) The calculated band structure of Ca3PbO along the 
high symmetry lines in the Brillouin zone. The location of one of 
the Dirac points is marked with an arrow. 



basis of the obtained results. 

The crystal structure of CaaPbO is shown in Fig. 1(a). 
It belongs to the space group PmSm 11 ' and is an inverse 
perovskite that possesses an O atom surrounded octahe- 
drally by Ca atoms. Figure 1(b) shows the band struc- 
ture of CaaPbO obtained within the first-principles cal- 
culation using the WIEN2k package, 12 ' in which the full- 
potential augmented-plane-wave method is implemented. 
The spin-orbit coupling is taken account of within the 
muffin-tin sphere of each atom via the second variational 
step. 13 ' The required crystal parameters are taken from 
experimental results. 11 ' 20 x 20 x 20 fc-points, which re- 
sult in 220 fc-points in the reduced Brillouin zone, are 
employed in the self-consistent cycle of our calculations. 
The charge density obtained with this number of fc-points 
is used to calculate band structures on the finer momen- 
tum meshes required to confirm the existence of Dirac 
electrons. An overview of the calculated band structure 
shown in Fig. 1(b) is as follows. First, the band at ap- 
proximately — 8 eV and the bands between —4 and eV 
mainly originate from Pb 6s and Pb 6p orbitals, respec- 
tively. On the other hand, the bands between —6 and 
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Fig. 2. (Color online) Dispersion relations around one of the 
Dirac points for (k x , k y , k z ) where 0.04 < k x < 0.2, —0.08 < k v < 
0.08, and k z has several values. The wave numbers are in the unit 
of 27r/ao, with ag being the lattice constant. 




Fig. 3. (Color online) (a) Naming convention of the three Ca 
sites and schematic view of the d-orbitals included in the basis set. 
(b) Definitions of hopping integrals. In both (a) and (b), the Pb 
sites are at the corners of the cube and the Ca sites are at the 
center of the surfaces. 



—4 eV originate from O 2p orbitals, and the highly en- 
tangled bands between and 6 eV mostly originate from 
Ca 3d orbitals. These observations mean that Pb 6s and 
O 2p orbitals are completely filled, and Pb 6p orbitals 
are almost completely filled, which is roughly consistent 
with Ca^+Pb^O 2 -. 

A Dirac point is a fc-point in the Brillouin zone around 
which Dirac electrons exist. We find that a Dirac point 
exists on the high-symmetry line connecting the T- and 
X-points, and its location is marked with an arrow in 
Fig. 1(b). Note that the marked point is exactly at the 
Fermi energy. Owing to the cubic symmetry, six Dirac 
points exist at the symmetrically equivalent points in 
the Brillouin zone. In order to clarify the dispersion rela- 
tions, three-dimensional plots of the band structure near 
the Dirac point are shown in Fig. 2. We can see a cone- 
like dispersion relation on the k x -k y plane when k z = 
[Fig. 2(a)] , and the gap between the valence band and the 
conduction band increases away from k z = [Figs. 2(b) 
and 2(c)]. These findings confirm the existence of three- 
dimensional Dirac electrons in the band structure. Ac- 
tually, even for k z — [Fig. 2(a)], the first-principles 
calculation gives a small gap of about 14 meV. Namely, 
the discovered Dirac electrons are massive. However, the 
estimated mass gap of about 14 meV is sufficiently small 
to observe physical phenomena associated with the linear 
dispersion of a Dirac electron. In fact, Fig. 2(a) shows 
that the linear (or sublinear) dispersion extends up to 
an energy higher than 0.2 eV, which is one order greater 
than the gap. 

In order to understand the emergence of the Dirac 
electrons in this material, we construct a tight-binding 
model describing the low-energy band structure. By an- 
alyzing the orbital-weight distribution in the obtained 
band dispersions 14 ) in Fig. 1(b), we find that the bands 
forming the Dirac electrons mainly originate from the Pb 
6p orbital and one type of Ca 3d orbital. It will become 
important later that an overlap exists between these p- 
and d-bands, i.e., the top of the p-bands (located at the 
T-point) has an energy of ~ 0.3 eV while the bottom 
of the d-bands (located at the T-point) has an energy 
of ~ —0.1 eV. Among the Ca 3d orbitals, those relevant 
to the Dirac electrons are the d y 2_ z 2 orbital on the Cal 
site, the d z 2_ y 2 orbital on the Ca2 site, and the d x 2_ y 2 
orbital on the Ca3 site, all of which are equivalent in the 
cubic symmetry. Figure 3(a) shows the definition of the 
Cal, Ca2, and Ca3 sites and the description of the above 



three orbitals. Based on these results, we use the 12 local 
orbitals 

\Px<j), \p y o), \p z cr), \dio), \d 2 cr), \d 3 a), (1) 

(a =f, I) as a basis set for our tight-binding model, 
where di, d 2 , and d 3 represent the above-mentioned d- 
orbitals on the Cal, Ca2, and Ca3 sites, respectively. As- 
suming the nearest-neighbor and next-nearest-neighbor 
hoppings shown in Fig. 3(b), we have the tight-binding 
Hamiltonian 

H = X] X e V C lpa^ V a + X e dcl dr7 C rda 
(7 rp a rd 

+ E E C^-r04«*<V«'* + 2A (2) 

o rar'a' r 

where indices a and a' represent p x ,y,z or di,2.3, and the 
fourth term represents the spin-orbit coupling for Pb 
atoms. Owing to the Pbp-orbital character, \p x f), \p y t)> 
and \p z l) are mixed. Here the spin-orbit couplings for 
Ca atoms are neglected for simplicity. The Fourier trans- 
formation of eq. (2) gives 

£ = EE f ~*'( fe ) c L<w. (3) 

k aa' 

After carrying out the Fourier transformation, we trans- 
form the matrix elements by attaching the momentum- 
dependent phase factors to the basis orbitals as 
\p x , y , z cr) -> e^+ k y+ k ^/ 2 \p XiViZ a), \d,a) -t e ik */ 2 \dxa), 
\d 2 a) -> e lk y/ 2 \d 2 a), and \d 3 a) -> e lk ^ 2 \d 3 a). This 
transformation makes the expressions for the matrix el- 
ements simple, and the transformed basis and matrix 
elements are used in the following. 

Figure 4(a) shows the band structure of our tight- 
binding model obtained using the parameter set e p — 
-1.25, t p i = 0.20, t p2 = 0.12, t p3 = 0.06, e d = 2.15, 
t d0 = 0.48, t dl = 0.22, t d2 = -0.22, t dp = 0.22, and 
A = 0.4 (in units of eV). Although there are some differ- 
ences between the first-principles bands and the tight- 
binding bands, the latter bands reproduce the impor- 
tant features (including the Dirac electrons) of the former 
bands having large Cal d y 2_ z 2 y Ca2 d z 2_ y 2, Ca3 d x 2_ y 2, 
and Pb p-orbital weights. In this simplified model, how- 
ever, there is no gap at the Dirac point, i.e., we obtain 
massless Dirac electrons. It is necessary to include Ca d xy 
or dzx/yz orbitals to have a finite mass term. (Details will 
be published elsewhere.) 
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Fig. 4. (Color online) Band structures obtained from the tight- 
binding model [eq. (2)] with (a) and without (b) the hybridization 
between p-orbitals and ci-orbitals. 



In order to obtain a clearer view on the emergence 
of the Dirac electrons, we derive a low-energy effective 
Hamiltonian of the above tight-binding model and show 
that it really results in a Dirac Hamiltonian. Firstly, by 
introducing a new basis set, we eliminate the uppermost 
and lowermost bands, which are well separated from the 
bands forming Dirac electrons. Since the Dirac point is 
close to the T-point, we use the new basis set, which con- 
sists of the eigenstates at the T-point. At the T-point, 
we can treat the p- and d-orbitals separately since the 
hybridization between them becomes zero. For the p- 
orbitals, the eigenvalues are e p o + A and e p0 — 2A, where 



e p0 



2t p i + 4i„2 + 4£ D 3 . We find that the former cor- 



<<p2 - 



responds to the states forming the Dirac electrons and 
the latter to the lowermost band. The energy splitting 
between them is caused by the spin-orbit coupling, and 



the former has a total angular momentum of j 
The new basis set for j — 3/2 states is 

|p f |) = -(b, t)+i\Py t»/V2, 



\p 
\p 
{Pi 



3 1 

2 2 



= -{\ Px i)+ l\ P y i)-2\ Pz t))/V6, 
= (\ Px t) - i\ P y t) + 2| P , I)) /TO, 

= (\ PX i)-i\ P yi))/V2. 



3/2. 

(4a) 
(4b) 
(4c) 
(4d) 



For the d-orbitals, the eigenvalues are £do~ 4£do and £do + 
8£do, where edo = (-d + 2tdi+^td2- The former corresponds 
to the bands near the Fermi energy and latter to the 
uppermost band. Explicitly, the wave functions for the 
former are 



\d[a) 



\d' 2 a) 



= (\d 1 a)-\d 2 a))/V2, (5a) 

= (|di(7) + |dao-)-2|d3cr))/V6. (5b) 

Performing a unitary transformation from the basis set of 
eq. (1) into that containing eqs. (4) and (5), we transform 
the 12 x 12 matrix, in eq. (3) into £' k . Then, we keep 
only the matrix elements between the above eight states 
in eqs. (4) and (5), reducing the matrix £' k into an 8 x 8 
matrix. Finally, concentrating on the Dirac point on the 
fc z -axis, we expand the matrix elements with respect to 
k x and k v up to the first order in k x and k y . After some 



algebra, we obtain 
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where k± = k x 



ci/y/3, c' 2 = C2/V3, and 



ci = it dp cos(fc 2 /2), c 2 = it dp {2 + cos(k z /2)) /V3, 



(6) 



(7) 



c 3 = — 8^d P sin(fc z /2)/\/3. 
Note that now the new basis set is {\p3 3), \p3 1), |p 3 1), 
|P||), K t), 14 t>, \d[ I), 14 I)}- Note alsV that 
and gf d depend on k z . 

In the limit of k x — k y — 0, all of the off-diagonal 
elements in eq. (6) vanish except for the elements rep- 
resented as C3. The finite C3 induces the strong band re- 
pulsion between bands from and \d[ f) (second 

and fifth rows), and also between bands from \p3i) and 

2 2 

\d[ l) (third and seventh rows). As a result, the bands 
originating from these states are pushed away from the 
Fermi energy, and thus we can neglect these states, keep- 
ing only |^»3 5 ) , 14 t)j an d 14 i) m the follow- 
ing. Then, suppose that g pp — g 2 d is satisfied at some 
k Z Q. When g\ p = g dd holds, we can expand g pp and 
g dd as g pp = —c p Sk z + eo and g dd = c d Sk z + e , where 
6k z =k z - k z0 and e = gf (= g dd ) at k z = k z0 . Using 
these relations, the Hamiltonian can be written as 



£'k = {eo+ScSk z )l- 



f-c5k z 



V 





-cSk. 



c 2 kj, 


c5k z 



\ 

-c 2 fc_ 



cSk z ) 



(8) 



with c = (cd + c p )/2 and 8c = (cd — c p)/2, which gives 
a tilted massless Dirac Hamiltonian. The Dirac point is 
(0,0, k z o) in this case. 

A very simple and intuitive reason for the emergence 
of the Dirac electrons is obtained from the above analy- 
sis. Figure 4(b) shows the energy dispersion of the tight- 
binding model, where td P is artificially set to zero. This 
corresponds to the case that there is no hybridization be- 
tween the p- and <i-orbitals. In this case, the p-bands (ri- 
bands) have a usual hole (electron) Fermi surface around 
the T-point. Although the p- and e?-bands overlap with 
each other, this only leads to a usual band crossing. If we 
turn on the hybridization td p , it generally causes a band 
repulsion between the p- and d-bands. However, owing 
to the symmetry of the involved orbitals, this band re- 
pulsion does not act on the T-X line. As a result, the 
band crossing remains at an isolated point in the Bril- 
louin zone, leading to the emergence of a Dirac electron. 
Note that the existence of an overlap between the p- and 
G?-bands supports the assumption g pp = g 2 d in the previ- 
ous paragraph. Thus, we can state that the overlapping 
between the p- and d-bands is essential for Dirac elec- 
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Fig. 5. (Color online) Allowed patterns of the phase of the wave 
function for the d-orbitals when k x = k y = 0, seen from the pr- 
axis direction. The dark and light robes represent the positive and 
negative values of the wave function, respectively. 



trons to appear in this model. 

Next, we discuss the reason why the band repulsion 
does not act on the T-X line. On the fc z -axis, which is 
equivalent to the T-X line, the ci-orbital wave function 
can be classified as shown in Figs. 5(a)-5(c). Actually, 
\d[a) corresponds to the wave function in Fig. 5(c), while 
\d' 2 (j) is a superposition of those in Figs. 5(a) and 5(b). 
From these figures, we can easily see that all of the hy- 
bridization between the Pb p-orbitals and the wave func- 
tions in Fig. 5 vanish except for that between the p z 
orbital and the wave function in Fig. 5(c). This is the 
reason why we have a finite off-diagonal matrix element 
C3 for | die), while the off-diagonal matrix elements for 
\d' 2 o) vanish as \k± \ — > 0. The vanishing matrix elements 
imply that the band repulsion does not act. 

Note that, although the above graphical argument is 
heuristic and intuitive, the spin-orbit coupling is not 
taken into account in this argument. By analyzing the 
irreducible representation of the bands obtained in the 
first-principles calculation, 14 ) we find that the two bands 
forming the Dirac electrons fall into the same irreducible 
representation on the T-X line. As a result, there is a 
small hybridization between them, leading to a small 
mass (gap) term. The arguments based on Figs. 5(a)-5(c) 
give a reason for the smallness of the induced mass gap. If 
the spin-orbit coupling is neglected, the two bands form- 
ing the Dirac electrons belong to different irreducible rep- 
resentations on the T-X line, so that a gap does not open. 
Namely, the spin-orbit coupling is necessary to make the 
mass term finite, even if d zx / yz or d xy orbitals are in- 
cluded in the tight-binding model. 

We discuss some points hereafter. The first is about 
related materials. We studied the band structures of 
the family of CaaPbO, namely, Sr 3 PbO, Ba3PbO, and 
CaaSnO. (Figures are not shown here.) In the calcula- 
tion, the lattice constants of these materials are again 
taken from experimental results. 11 ' We find that all these 
systems have similar band structures and Dirac elec- 
trons. We point out here that the overlap of the p- and 
(i-bands, which is essential for the existence of Dirac elec- 
trons, becomes larger in the order CasPbO — > SraPbO — > 
Ba3PbO. However, in Ba3PbO, some other (non-Dirac) 
bands appear at the Fermi energy, which may mask the 
properties of the Dirac electrons. For CaaSnO, the over- 
lap of the p- and <i-bands becomes smaller than that in 
CaaPbO, but it still exists and gives Dirac electrons at 
the Fermi energy. Since the spin-orbit coupling is smaller 
for Sn than for Pb, the estimated mass gap is actually as 
small as 4 meV, compared with 14 meV for CaaPbO. 



From this result, we propose that the series of alloys 
Ca3(Pbi_ x Sn a; )0 will provide a method of controlling 
the mass term of Dirac electrons. 

Compared with Bi, which has three-dimensional mas- 
sive Dirac electrons in its low energy band structure, the 
present system has a simple structure whose Fermi en- 
ergy is located in the gap of the Dirac electrons. In Bi, 
on the other hand, there is a complication due to its 
semimetallic properties. 15 ' As a result, fine-tuning of the 
doping or pressure is necessary for studying the Dirac 
electrons in Bi. Ca3PbO and its family will be inter- 
esting and important materials for exploring the three- 
dimensional Dirac electrons. 

Responses to applied magnetic fields will also be in- 
teresting. Neglecting the small gap, a Dirac electron has 
a linear dispersion, leading to a peculiar Landau level 
structure. 16 ' Another possible interesting response is a 
large orbital diamagnetism, which is observed in Bi. 17 ' 18 ' 
The orbital diamagnetism takes a maximum when the 
Fermi energy is inside the gap, namely, a large electronic 
response is observed even though there are no Fermi sur- 
faces. The interesting interband effects in the quantum 
transport phenomenon are being extensively studied. 19 ' 

If a small number of electrons (or holes) are doped to 
this material, it becomes a metal with six tiny Fermi sur- 
faces, or six valleys. It will be interesting to consider the 
possibility of a three-dimensional analogue of valleytron- 
ics, which has recently been developed for graphene. 20 ' 21 ' 
It is also interesting to consider the possibility of sponta- 
neous valley symmetry breaking, in which the occupancy 
of the valleys spontaneously becomes imbalanced. 

In summary, we have performed the band structure 
calculation of CaaPbO and its family, and found that 
three-dimensional massive Dirac electrons exist at the 
Fermi energy. A tight-binding model that describes the 
essence of the emergence of Dirac electrons was success- 
fully constructed. By analyzing this model, it is con- 
cluded that the symmetry of the crystal and the involved 
orbitals play important roles in sustaining Dirac elec- 
trons in the band structure. 
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